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1 Introduction 

In this paper we describe the geodesies of a left-invariant sub-Riemannian 
metric on the three-dimensional solvable Lie group Sol. 

This group is widely known in geometry, because it allows compact 
quotient-spaces and it gives one of the Thurston three-dimensional geome- 
tries pQ. By the classification theorem of Agrachev-Barilari [2 J there are 
invariant sub-Riemannian geometries realized on four solvable non-nilpotent 
Lie groups: SE(2), SH(2), SOLV~, and SOLV+. 

In this classification, our geometry corresponds to the case SOLV~: 

Sol = SOLV~. 

The case of SOLV + we shall consider separately. 

Various aspects of the integration of geodesic flows on sub-Riemannian 
manifolds has been widely studied (see, for example, [31 HI [5j Q2]). Note, 
that the geodesies of other three-dimensional nonsolvable or nilpotent sub- 
Riemannian geometries have been described recently in terms of elementary 
functions [TIE]- In our situation it is necessary to use elliptic functions. 

We thank LA. Taimanov for posing the problem and Ya.V. Bazaikin for 
helpful discussions. 
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2 Basic definitions 



2.1 Geodesies of sub-Riemannian manifolds 

Let M n be a smooth n-dimensional manifold. A smooth family of /c-dimensional 
subspaces in the tangent spaces to points of M n 

A = {A(q): A(q) G T q M n V?eM n , dim A(g) = A;} 

is called completely nonintegrable, if the vector fields tangent to A, and all 
their iterated commutators generate the tangent bundle TM n : 

span { [ft, [. . . [fm-ljm] .-.]](?): fM G A (p) Vp £ M n , m = 1, . . . } = T q M n . 

Sometimes this distribution is called completely nonholonomic. 

A two-dimensional distributions in a three-dimensional manifold is com- 
pletely nonholonomic if and only if 

span{/i( g ),/ 2 ( 9 ), [fi(q),f 2 (q)]} = T q M 3 , 

where at every point q the vectors fi(q) and /2(g) form a basis in A(q). 

Let gij be a complete Riemannian metric on M n . A triple (M n , A, t^,-) 
is called a sub-Riemannian manifold. A Lipschitz continuous curve 7 : 
[0,T] -> M n is called admissible if j(t) G A( 7 (t)) for almost all t G [0,T]. 
The length of this curve is equal to 

The distance between two points on the manifold is defined by formula 

d(qo,qi)= inf 

7esz 90 , 9l 

where ^ 90l91 is the set of all admissible curves connecting points go an d 
q\. This function d(-, ■) is called the sub-Riemannian metric on M n . A 
geodesic of this metric is an admissible curve 7 : [0, T] —> M n , which locally 
minimizes the length functional £(7). 

Geodesies of sub-Riemannian metrics satisfy the Pontryagin maximum 
principle (see, for instance, [5]), which we formulate below. Let /1, •••,/& 
be vector fields which are tangent to A and span A at every point of M n 
(or of a domain of M n ). 

The Pontryagin maximum principle is as follows: 
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• Let M n be a smooth n-dimensional manifold. Let us consider for 
Lipschitz continuous curves the following minimum problem 

k ,.T k 

q = y ^2uifi{q), Ui£R, / VNf(i)di — > min, q(0) = q , q(T) = q 1 
i=i Jo i=i 

with a fixed T. Let us consider the mapping H : T*M n xlxtf^l, 
given by the function 

k k 
H(q,\,p ,u) := {X,^2uifi(q)} + Po^u 2 i . 

i=i i=l 

Jf a curve q(-) : [0,T] -)• M n with a control u(-) : [0, T] ->■ M fc is 
optimal, then there exists Lipshitzian covector function A(-) : t G 
[0,T] ^ A(t) G T* (t) M n , (A(t),p ) / and a constant p < such 
that 

*J ?(*) = ^-(?(t),A(t),po, «(*)), 
dT-L 

n) A(t) = — Q-(q(t),\(t),Po,u(t)), 
Hi) —(q(t),X(t),p ,u(t)) = 0. 

A curve q(-) : [0, T] —> M n , satisfying the Pontryagin maximum principle 
is called an extremal (curve). To such a curve there corresponds a set of 
pairs (\(-),po). The type (normal or abnormal) of an extremal depends on 
the value of po: 

• if po 7^ 0) then the extremal is called normal; 

• if po = 0, then the extremal is called abnormal; 

• extremal is called strictly abnormal if it is not projected (on M n ) onto 
a normal extremal. 

For a normal extremal we may put po = —\. 

Normal extremals are geodesies [5J. In the contact case, when at every 
point the distribution A coincides with the annihilator of the contact form 
on M n , there are no nontrivial abnormal extremals (this fact is indicated 
in [7j). In the case, when the space of vector fields on a manifold is gen- 
erated by vector fields tangent to the nonholonomic distribution and their 



3 



commutators, there are no strictly abnormal extremals [5]. We recall, that 
an extremal is called strictly abnormal, if its projection on M n does not 
coincide with the projection of any normal extremal. Both of the above 
statements apply to three-dimensional sub-Riemannian manifold M 3 . 

By hi), Uj = (X(t), fi(t)) and a curve q(-) : [0, T] — > M n is geodesic if and 
only if it is the projection onto M n of a solution (X(t),q(t)) of a Hamiltonian 
system on T*M n with the following Hamiltonian function: 



€ M n , AG T*M n . 



The Hamiltonian H is constant along any solution of the Hamiltonian sys- 
tem. Moreover, H = ^ if and only if the geodesic is length parameterized. 

2.2 Elliptic functions. Jacobi functions 

Let us recall some necessary facts of Jacobi elliptic functions. The integrals 

dx 



f 

Jo 



/o VU-^X 1 -*^ 2 )' 
and 

x Vl - k 2 x 2 , 
— dx, 

are called elliptic integrals of first and second kind, respectively, in the nor- 
mal Legendre form (see |10l I13j). where k (0 < k < 1) is the modulus of 
these integrals, k! = Vl — k 2 is the additional modulus. By the substitution 
x = sin ip these integrals reduce to the normal trigonometric form 



da f sm ^ dx 



V(l - k 2 sin 2 a) Jo ^(1 - x 2 )(l - k 2 x 2 ) ' 



ff I raimp Vl — k 2 X 2 

E(<^, k) = v 1 — k 2 sin 2 a da = — dx 

Jo Jo a/1 — ^ 



vr 



X 



Consider an integral of first kind in the normal trigonometric form 

<"f d(p 

v 



yi — k 2 sin 2 (p 
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Now consider the upper limit as a function of v. This function is denoted 
by 

<p = am (v, k) = am v 

and is called the amplitude, and this process is called an inversion of the 
integral. Thus, next functions: 

sin ip = sin(amu) = snu, 

cos (p = cos(amu) = cnu, 



Aam v = y 1 — k 2 sin 2 if = \/l — k 2 sn 2 v = dn v 
are called the Jacobi functions and are related by 



sn v + cn v = 1, dn v + k sn v = 1. 



By derivation, we obtain 

dsnv 



dv 

dcnv 

dv 
ddnv 



dv 

and conclude, that 

2 



cnudnt;, 

— snvdnw, 
-k 2 sn vcnv 



I dsn v \ , 2 \ /i j 2 2 , 

1 — sn v)(l — k sn v j, 



\ 2 

= (1- cn 2 v)(k' 2 + k 2 cn 2 v), (3) 



<icnt; x 



dv 

\ 2 

/2^ 



/ ddnv^ 



V (if 



(l-dn^)(dn^-^). 



The first equation of ([3]) implies that sn v is the inversion of the first kind 
elliptic integral in the normal Legendre form 

dX (4) 



o V(l - x 2 )(l - k 2 X 2 ) 
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Prom the second and third equations obtained that cn v and dn v are the 
result of conversion of next functions 

f cnv dx 
Ji ^/{l-x 2 ){k 12 + k 2 X 2 ) 

r , ix <6) 

Jl y/(l ~ X 2 ){ X 2 - W) 

All Jacobi functions are periodic. Note that the function snu is odd, 
but cnu and dnu are even, therefore we assume, what in the two last in- 
tegral, when the functions cnt) and dnv pass through the critical points, 
respectively changes a sign of radical. 

3 Sub-Riemannian problem on the group SOLV 

Let us consider the three-dimensional Lie group SOLV - formed by all ma- 
trices of the form 

e~ z x \ 
e z y ,x,y,zeR. 
1/ 

Its Lie algebra is spanned by the vectors 






ei = U U U , e 2 = U U l , e 3 



meeting the following commutation relations: 

[ei,e 2 ] = 0; [ei,e 3 ]=e 1 ; [e 2 ,e 3 ] = -e 2 . 

We take a new basis 

ai = ei + e 2 ; a 2 = e\ - e 2 ; 03 = es, (7) 

in which the commutation relations take the form 

[ Ol ,a 2 ] = 0, [ai,a 3 ]=a 2 , [a 2 ,a 3 ]=ai. 

Let us consider a left-invariant metric on SOLV - , which is defined by its 
values in the unit of the group: 



(ei,ej) = 5,, 
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The Lie group SOLV~ is diffeomorphic to the space M 3 . Indeed, x,y,z 
are the global coordinates on SOLV - and they also may be considered as 
global coordinates on M 3 . The tangent space at each point of SOLV - is 
spanned by matrices of the form 



& = I | , dy 



which are the left translations of the basic vectors: 

L q *{e\) = e~ z d x , L qSf (e 2 ) = e z d y , £ g *(e 3 ) = d z . 
Since the metric is left-invariant, we have 







' -e~ z 








d z = ( 





e 2 







. o 









e 2z 











e -2z 
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9ij{x,y,z) = 

For the basis oi, 03 we have 
L q *(ai) = e' z d x + e z d y , L g *(a 2 ) = e~ z d x - e z d y , L q *(a 3 ) = d z 
The inner product takes the form 




(L q *(ai),L q *(aj)) = (a*, a,) = 2 (8) 



In this paper we study the sub-Riemannian problem on the three-dimensional 
Lie group SOLV - defined by the distribution A = span{ai,a3} with metric 

©• 

Let G = SOLV , Q be its Lie algebra with the basic vectors 01,02,03 
([7]). We split the Lie algebra Q into the sum sum pQ)k, where p = 
span{oi,03}, A; = span{a2}. 

Let us consider a two-dimensional left-invariant distribution A = spanjai, 03} 
in TG, and a left-invariant Riemannian metric (|8|) for which the spaces p 
and k are orthogonal, i.e., the metric tensor splits as follows: 

9 = (9ij) = 9 P + 9k- 
Let us introduce a parameter r and consider the metrics 

9r = 9 P + rg k - 
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Every such a metric together with A defines the same sub-Riemannian man- 
ifold because only the restriction of the metric onto A is important. 

However the Hamiltonian function for the geodesic flows of these metrics 
depends on r: 



H(x,p,t) = -g^{x)piPj, 



where g^g^ k = Sf. We have 



9r,ij 



( ^e 2z V \ 



V o 







1 / 



9? 



l-r 1±t c 2z q 



2r 







1 J 



The Hamiltonian function H for the normal geodesic flow of the sub-Riemannian 
metric is obtained from H(x,p,r) in the limit 



t — y oo, 



and we derive 



1 1 

H(x,y,z,p x ,p y ,p z ) = -e~ 2z pl + -p x p y + j e lz p l y + ^pj. 



1 2z 2 , 1 2 



(9) 



The Hamiltonian equations = ^ , pj = — ^ take the form 



an 



■ 1 -9, 1 

x = -e zz p x + -pj,, ^ = 0, 



■ 1 9, 1 

V= ^e zz p y + -p x , 



z=p z 



Py = 0, 
1 

Pz 



-2zJ2 



a 2zJ2 



(10) 



These differential equations can be derived from the Pontryagin maxi- 
mum principle. The corresponding Hamiltonian takes the form 



1 



H(x,y,z,p x ,p y ,p z ,p ,u 1 ,u 3 ) = -^={uip x e z + uip y e z ) +u 3 p z +p {u\+ul), 



where po = — -, u\,u 3 are control functions. 
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The system (|10p has three first integrals: 



h = H, h=p x , h=P y , 

which are functionally independent almost everywhere, and therefore the 
system is completely integrable. 

Since the flow is left-invariant as well as the distribution A and the 
metric, without loss of generality we assume, that all geodesies originate at 
the unit of group, that is, we have the following initial conditions for the 
system (fTU|) : 

x(0) = 0, y(0) = 0, z(0)=0. (11) 



In the sequel, we put 



u 1 Px Py , 

H = -, — = = a, —7= = b. 



2' V2 ' y/2 
By substituting these expressions into ([9]), we obtain 

1 = ( e -z a + e * b f +p l (12) 

which implies 

p z = ±\J 1 — (e _z a + e z b) 2 . 

By substituting this expression to the third equation of (jlOp we obtain 
equation for the temporal variable t for positive values of p z 

" Z (13) 



1 - (e~ z a + e z bf 



If p z < 0, then all calculations will be similar, but with contrary sign. 
Let us make the change of variables 



u 

and rewrite (fT3l) as 



(14) 



u 2 — (a + bu 2 Y 



The last expression is not integrated in terms of elementary functions and 
defines an elliptic integral, except of special cases, when this elliptic integral 
degenerates. These cases will be discussed below. 
Let us consider the generic case o/O and 6^0. 
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The subradical expression in (|14p haves discriminant D = 1 — Aab > 
accordingly (fT2j) . 

D = if and only if p 2 = accordingly system (fT0|) and equation (fT2|) . 
That case is degenerative. 

Thus, if D > < — V then there exist o\ and cr|, such that the 

following 

u 2 - (a + bu 2 ) 2 = -b 2 u 4 + (1 - 2ab)u 2 - a 2 = 



= -b\u 2 - a 2 )(u 2 - a 2 ) = a\b 2 U-t\(± °A 

V °i/ V°i °i 

and 

''1,2 " ^2 



2 1 - 2a& ± VI - 4a6 

*1,2 = ST5 ( 15 ) 



Put 

u 

io = — (16) 



and rewrite (|14j) in the following form 

1 

o\b 



;, f " ; , 

(7? 



We apply the Jacobi elliptic function (jBJ) in order to inverse this integral: 

cdn(<7iW) dw 



o\bt 



2 

where k' 2 = -§■. Therefore 



10 = dn (a\bt, k), 

where 

fc 2 = l-4- (18) 
CT i 

By inverting f)16[) . putting u = e z , and keeping in mind the initial condition 
(fTT|) and the equality dn (0, k) = 1, we obtain 

z(i) = lndn (a\bt, k). 
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By substituting this expression into the first equation of (jlOp and integrating 
it in elliptic functions (seepjj]), we derive: 



x(t) = 



/c 2 sn (<ti 6i) cn (a\bt) 1 . , . ,.\ 

' +- n5 E(am(<7 1 W), fc) +bt 



+ C, 



where E(x,k) is elliptic integral of second kind 

Since sn (0, k) = 0, cn (0, k) = 1, am (0, fc) = and E(0, fc) = 0, we have 
C = 0. From the second equation of this system we conclude that 

y(t) = —= I — E(am(<Ti6t), k) + at J + Q, 
V2 \o-l / 

with Q = const. By (fTT]) . we compute that Q = 0. 

Let us now consider the cases, when the elliptic integral (j!4h degenerates: 

1. a = 0, 6 = 

2. a = 0, 6/0 

3. 6 = 0, o^0 



4. 5 = ol) 



1 



We consider its successively: 

1) a = 0, 6 = 0. From the equations (fT0"|) . (fTTj) and (fl~3|) is cleaner, that 



x(t)=0, y(t) = 0, z(t)=t. 

2) a = 0, 6 / 0. We have p x = 0, p y = y^fc. 
The equation (jl4[) is rewritten as 



(19) 



6 2 u 4 



By integration and transformation the resulting expression by the inverse 
change of variable, we obtain u = e z : 

Cbe z 



l + Vl-b 2 e 2z ' 



where C = const and C > 0. The last expression together with the initial 
condition (fTTj) implies 



C 



li 



and we derive that 

2Ce l 



6(C 2 + e 2 *)' 

i.e., 

z{t) = ln b{c 2 + e*y 

which after substituting the formula for C takes the form 



,< t )-ln. 2 d + ^V 



2(1 + VT^W) -b 2 + b 2 e 2t 
By the first equation of (110p . we have 

b 

x(t) = - / =t, 



and the second equation of (flOj) together with (|20j) and (fTTj) implies 



6(C 2 + e 2i ) 6(C 2 + 1)' 
Finally in the case 2) we have the explicit formulas for solutions: 
b 

x(t) = —t 
V2 



V2 (2(1 + v 7 !^ 2 ) - 6 2 ) \/2 (2(1 + Vl^h 2 ) - b 2 



y(t) = r, — . ,/o, + 

z(t) = In 



26(1 + VT^¥) - 6 3 + b 3 e 2t 26(1 + VT^¥) 

2(l + v / T^6 2 ) e* 



2/ 



2(1 + v 7 !^ 2 ) - 6 2 + 6 2 e 

3) 6 = 0, a^O. We have p y = 0, p x = \/2a. 
The equation (|14p takes the form 



— or 



Let us put « = e z and derive 



e * = ( e 2 + Ve 22 - a 2 ) C, 
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where C = const and C > 0. The last expression together with (jlip implies 

C= . (22) 

from which we obtain 

where C is given by (122j) . As in the case 1) we derive from the first two 
equations of (fTUj) that 



x/2aC 2 V2aC 2 a 

l(i) = -^cv + ITcv' (23) 



Finally we obtain 



= , ^ — , — + ^' 



e 2t 



2(1 + VT^a 1 ) - a 2 + a 2 2(1 + \/T 



a 2 ) 



y(t) = ^*' (24) 



zt = In + ^ i— 

W \ v 2(l + v / T^)e* 2 

4) L> = (aft = 1 
Note that, the formula (1121) by (jllh is rewritten as 

(a + 6) 2 +p 2 = l, (25) 

which means that 

|o + 6|<l. (26) 

Then it is clear that a = b = ^ or a = b = and for these values the 
equation ([25]) implies that p z = 0. Therefore solutions to (fTOj) in the case 4) 
are linear: 

x(i) = — =, y(t) = —=, z(t) =0, a = b = -; 

i t i \ ZI ) 

x(t) = --j=, y(t) = --j=, z(t) = 0, a = & = --. 

Thus we have the following 
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Theorem 1 In a generic case the normal geodesies (with the initial condi- 
tion Ul\)) are described by the formulas (for p z > 0): 



x(t) 



ak 2 sn (a\bt) en {a\bt) 



+ 



oE(am (o~\bt), k) 




b 



V2o-!bk> 2 y/l-k 2 sn 2 (<ribt) 



V2a!bk' 2 



y(t) 



E(am (o~ibt), k) 



+ 



a 

V2 



(28) 



z(t) 



In dn {o~\bt) 



where the parameters o~\ and k are determined by a and b (ab < 4) via U5\) 
and pap. 



In the degenerated cases 1~4 the normal geodesies (with the initial con- 
dition hll}) ) are described in terms of elementary functions by the formulas 



Notice that normal geodesies in the theorem are parameterized by a, b. 
Constants k, o\ are defined by a, b as explained before. 

The qualitative behavior of generic normal geodesic is quite complicated. 
Figures 1 and 2 show parts of the geodesic spheres of radius of 0.15 and 0.25 
(a scale on each figure itself; axis z is exponentially scaled). A grid on the 
spheres corresponds to two parameters 9 and where 6 is angle of the 
initial vector of the geodesic with respect to the axis x and [i is the initial 
acceleration value x + y along the geodesic, i.e. fi can be interpreted as the 
acceleration with which the geodesic is drawn out of the starting point. In 
Figures 6 varies from 7r/6 to 57r/6 (part of the sphere of — it/ 6 to — 5ir/6 is 
obtained as a mirror). Parameter jj, varies from —45 to 45. 

On this grid it can be seen only the qualitative behavior of a sphere 
with increasing radius. The drawings practically do not show those parts of 
spheres, which are too fast going to infinity, as well as those, which coincide 
to the geodesies, changing too quickly the direction. We can see that part 
of geodesies starting at small angle to plane x, y goes to large values of 
coordinates x, y very quickly, even for not large values of parameter fi. For 
sufficiently large 9 and mean values of \fi\ geodesies deviate not too much 
from plane x = y, but if \n\ increases the deviation from this plane begins. 
Assuming exponential scale of axis z we see that coordinate z increases much 
more slow than x and y. 



EW, M>, (M>> and & 
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Figure 2: r=0,25 
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